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Two-point functions of four-dimensional simplicial quantum gravity 

W. Beirl, H. Markum and J. Riedler* 

Institut fur Kernphysik, Technische Universitat Wien, A-1040 Vienna, Austria 

We investigate the interaction mechanism of pure quantum gravity in Regge discretization. We compute 
volume-volume and link-link correlation functions. In a preliminary analysis the forces turn out to be of Yukawa 
type, at least on our finite lattice being away from the continuum limit. 
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1. INTRODUCTION 

Regge calculus gives a description of Gen- 
eral Relativity using simplicial approximations 
of space-time geometries [1]. A Regge skeleton 
consisting of a set of four-simplices is completely 
specified by its incidence matrix and by the link 
lengths. Here we use skeletons with a fixed inci- 
dence matrix restricting the topology to a four- 
torus and consider the link lengths as the dy- 
namical variables of the theory. To define the 
Euclidean path integral one must choose a func- 
tional measure which ensures that the simplices 
satisfy the triangle inequalities and their higher 
dimensional analogs. 

The functional integral with the Regge- 
Einstein action is 
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where triangle areas A t , deficit angles b t , and 
four-simplex volumes V s are calculated from the 
squared link lengths q\. The function T is equal 
to one for Euclidean configurations and zero oth- 
erwise and the parameter a allows to consider dif- 
ferent types of measure. Scaling arguments sug- 
gest to set the cosmological constant A = a [2]. 

The phase structure of the theory (1) has al- 
ready been studied extensively [2-4]. We men- 
tion only that in spite of an unbounded gravita- 
tional action the system develops a well-defined 
phase with small negative average curvature for 
couplings [3 < f3 c . In order to facilitate numerical 
simulations a lower limit / = 10 -4 < <f> s has been 
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applied to the fatness of every four-simplex [4,5] 

V 2 



max ;es (gf) 



2. CORRELATION FUNCTIONS 
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To understand the interaction mechanism of a 
theory one has to study correlation functions. In 
case of the Regge approach the observables O(x) 
are associated with a vertex x and the corre- 
sponding expectation values 



G(d) = (O(x)O(y)) 



(3) 



should be measured at two points separated by 
the geodesic distance d = \x — y\ [2]. As ob- 
servables one can consider quantities like the lo- 
cal volume, the local curvature, the link lengths 
or the deficit angles. Although not all of these 
objects are diffeomorphism invariant in the con- 
tinuum their computation yields insight into the 
behavior of the Regge skeleton. 

For a first study we investigated connected 
volume- volume correlation functions 



Gv(d) = (V(0)V(d)) c 



(4) 



where V denotes the hypercubic volume element 
consisting of 24 four-simplices associated with 
each vertex. We performed MC simulations for 
several (3 values on lattices with 3 3 x 8 and 4 3 x 16 
vertices applying 100k - 300k and 30k - 70k 
sweeps, respectively. For asymmetric lattices it 
seems justified to use the index n labeling the ver- 
tices in the long direction instead of the geodesic 
distance d. This avoids the time consuming deter- 
mination of geodesic distances on the fluctuating 
simplicial lattice. 
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Figure 1. Volume correlations Gy(n) for several 
couplings (3 and A = a = 1 on a 4 3 x 16 lattice. 
The fits suggest a Yukawa like theory. 



Figure 2. Volume correlations Gv(n) for several 
couplings (3 and A = a = 0.1 on a 4 3 x 16 lat- 
tice. Screening masses decrease with decreasing 
cosmological constant. 



In Fig. 1 we present our results for the volume 
correlations in the case of a uniform measure and 
cosmological constant A = 1. We find a fast decay 
of the correlators for all investigated [3 values near 
the critical coupling f3 c > 0.108. In a preliminary 
analysis of the data we only tried an exponen- 
tial decay without any power law term. We men- 
tion that the corresponding mass lies in the range 
2 < m < 8 in inverse lattice units. The cosmolog- 
ical term is expected to lead to a non-vanishing 
screening mass, m ~ A [6]. Therefore, we low- 
ered A = a to 0.1 approaching a scale invariant 
measure. In Fig. 2 we display the volume corre- 
lators for several [3 values near f3 c > 0.14. Note 
the difference in the scale by a factor of 10. A fit 
to an exponential decay yields an effective mass 
2 < m < 4 being smaller than in the A = 1 case 
as naively expected. Since it was demonstrated 
in earlier work that one-point expectation values 
are practically identical in the range - a < 1 [4] 
it seems reasonable to compare directly the effect 
of the cosmological constant. 

To gain further information about the geom- 
etry of the simplicial lattice we computed the 
15 x 15 correlation matrix between the (squared) 



link lengths 

Qab(d) = (q a (0)qb(d)} c (5) 

and investigated its eigenvalue spectrum A a (d). 
The indices a, b label the different types of edges 
within a hypercube along the indexing scheme of 
Ref. [7]. Some care must be taken to interpret 
Qab(d) as a propagator. From perturbative cal- 
culations considering small fluctuations around 
flat space [7] it is known that the free inverse 
propagator on the lattice has zero eigenvalues (as 
well as in the continuum). One can argue that 
these zero modes have no extensions to a general 
curved lattice since diffeomorphism invariance is 
lost. The formulation of gauge transformations 
and gauge fixing are merely worked out for lat- 
tice gravity. Besides, the construction of physical 
two-point functions has been criticized even for 
continuum quantum gravity [8] but the fact that 
the lattice theory is formulated in a coordinate 
invariant manner at least addresses one of these 
objections [2]. 

In Fig. 3 we draw our results for the largest 
eigenvalue of the link correlations for uniform 
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Figure 3. Largest eigenvalue A(n) of the link- 
correlation matrix Q a i, for A = a = 1 on a 4 3 x 16 
lattice and the same couplings [3 as in Fig. 1. 



Figure 4. Largest eigenvalue A(n) of the link- 
correlation matrix Q a i, for A = cr = 0.1 on a 4 3 x 
16 lattice and the same couplings [3 as in Fig. 2. 



measure with A = 1. Because our data shows 
large fluctuations we simply fitted an exponential 
without power law behavior and got for the mass 
0.3 < m < 0.7. Fig. 4 gives similar link corre- 
lations for an approximately scale invariant mea- 
sure with small cosmological term A = a = 0.1. 
An exponential fit yields 1 < m < 3. 

3. CONCLUSION 

This work contains a first computation of 
two-point functions for four-dimensional simpli- 
cial quantum gravity. The preliminary results 
should be understood as a status report of a 
project considering different types of correlations. 
From the viewpoint of weak-field approximation 
with all the problems of a continuation to non- 
perturbative quantum gravity one might expect 
the curvature correlator as a very promising can- 
didate to carry the quantum numbers of the gravi- 
ton [2]. Our exploratory data analysis for the 
theory at finite lattice spacing indicates Yukawa 
forces with relatively large masses. The crucial 
question is the behavior of the mass spectrum in 
the continuum limit. 
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